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Abstract 

We analyze the second order QCD corrections to the fragmentation functions F^(x,Q 2 ) 
(k = T,L,A) which are measured in e + e~ annihilation From these fragmentation functions 
one can derive the integrated transverse (or), longitudinal (o"x) and asymmetric (a a) cross 
sections. The sum a tat = ot + ol corrected up to order c? s agrees with the well known 
result in the literature. It turns out that the order o? s corrections to the transverse and 
asymmetric quantities are small in contrast to our findings for F^(x,Q 2 ) and ol where 
they turn out to be large. Therefore in the latter case one gets a better agreement between 
the theoretical predictions and the data obtained from he LEP experiments. 
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1 Introduction 

Semi inclusive hadron production in e + e~ annihilation into a vector boson V (V = 7, Z) proceeds 
via the reaction (see H) 

e-(h, ai) + e + (Z 2 , a 2 ) - V(q) H(p, s) + "X" . (1) 

In the process above "X" denotes any inclusive final hadronic state and H represents either a 
specific charged outgoing hadron or a sum over all charged hadron species. The unpolarized 
differential cross section of the above process is given by 

d 2 cr H (x,Q 2 ) 3 2 ^ da$(x,Q 2 ) , 3 , 2 dof (x,Q 2 ) 

, , 7- = -(1 + cos 9) H7 sin 9 

ax a cos 9 8 ax 4 ax 

4 dx 

The transverse, longitudinal and asymmetric cross sections are given by a? , , and crjf re- 
spectively. The latter, which is due to parity violation, only shows up if the intermediate vector 
boson is given by the Z-boson and is absent when V — 7. The cross sections depend in addition 
to the CM energy Q also on the Bj0rken scaling variable x defined by 



S = 7?F' q 2 = Q 2 >0, 0<x<l. (3) 

Q 



In the centre of mass (CM) frame of the electron-positron pair this variable can be interpreted as 
the fraction of the beam energy carried away by the hadron H. The variable 9 denotes the angle 
of emission of particle H with respect to the electron beam direction in the CM frame. Before 
the advent of LEP1 the CM energies were so low (Q Mz) that a a could not be measured and 
no effort was made to separate ot from <jl so that only data for + were available. After 
LEP1 came into operation one was able to measure ot and ai separately 0, @. Moreover a a 
could be determined for the first time 0, M. The separation of or and ox is important because 
the latter cross section enables us to extract the strong coupling constant a s and allows us to 
determine the gluon fragmentation density D^(x, fi 2 )) with a much higher accuracy as could be 
done before. Furthermore the measurement of a a provides us with information on hadronization 
effects since the QCD corrections are very small. In the QCD improved parton model, which 
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describes the production of the parton and its subsequent fragmentation into a hadron H : the 
cross sections aj? (k = T, L, A) can be expressed as follows 



da£(x,Q 2 ) f 1 dz 



dx J x z 



vloliQ 2 ) {D»(*,n 2 )C s k!q (z,Q 2 /n 2 ) + Df(*,^ 
L z z 



<,(s,QV)} + E °?(Q 2 )D H NSJ (~ /) C^(z,Q 2 /^% (4) 



/=i 

for k = T,L. In the case of the asymmetric cross section we have 

7=i 



1 A n ? 

I — [ E A f (Q 2 ) D li(-> ^) c A S q («, Q 2 /V)] • (5) 



In the formulae above, which only hold for massless quarks, we have introduced the following no- 
tations. The functions Df(z, Li 2 ) (I = q, q, g), which depend on the factorization/renormalization 
scale (a, stand for the parton fragmentation densities. The singlet (S) and non-singlet (NS,A) 
combinations with respect to the flavour group SU(rif) are defined by 



D% s , q (z,Li 2 ) = D»(z^ 2 )+D^(z,Li 2 )-D»(z,Li 2 ), 
D% q (z,LC 2 ) = Df(z,LC 2 )-Df(z,LC 2 ). (6) 

The index q stands for the quark species and nf denotes the number of light flavours. The 
pointlike cross section crp and the asymmetry factor A' ' of the process e + + e~ — ► q + q can 

be found in ||, @. The total cross section, summed over all flavours is given by 0tot(Q 2 ) = 
J2qLi &q{Q 2 ) ■ The QCD corrections in eqs. © are described by the coefficient functions 
C r kl (k — T, L, A; I = q, g). Like the fragmentation densities they depend on the scale /i and they 
can be split into a singlet (r = S) and a non-singlet part (r = NS). From (||D we can derive the 
total hadronic cross section 



(?tot(Q 2 ) = 77 E / dx [ dcos6(x 
2 u Jo J-i v 



H u-l dxd cos 9 

y2\ 



with 



a T (Q 2 ) + a L (Q 2 ) , (7) 



cr k{Q) = -l^ / dxx 3 > fc = T,L 

2 x Jo dx 



where one has summed over all types of outgoing hadrons H. From the momentum conservation 
sum rule given by 

E/ dxxD?(x,Li 2 ) = 1, l = q,q,g, (9) 



2 



and Eqs. @, (§) one can derive 

<y k {Q 2 ) = ai°l(Q 2 ) f Q dxx[C s k>q (x, Q 2 //i 2 ) + l -Cl g {x, Q 2 /^ 2 )} . (10) 
Finally we also define the transverse, longitudinal and asymmetric fragmentation functions 

One observes that the above fragmentation functions are just the timelike analogues of the 
structure functions measured in deep inelastic electron-proton scattering. 

2 Order a 2 s corrected coefficient functions 

The coefficient functions C T kl corrected up to order a 2 receive contributions from the following 
parton subprocesses. In zeroth order we have the Born reaction 

(12) 

where "f (/ = q, q, g) denotes the detected parton which subsequently fragments into the hadron 
of species H. In next-to- leading order (NLO) one has to compute the one-loop virtual corrections 
to reaction (|T^) and the parton subprocesses 

V^"q"+q + g, (13) 

V^"g"+q + q. (14) 

After mass factorization of the collinear divergences which arise in the above processes one 
obtains the coefficient functions which are presented in Q]. The determination of the order 
a 2 contributions involves the computation of the two-loop corrections to (|i~2"D and the one-loop 
corrections to eqs. (|T3|) , (|14|) . Furthermore one has to calculate the following subprocesses 

V ->V +q + g + g, (15) 

V +q + q + g, (16) 

V ^"g" +q + q + q. (17) 

In reaction ([17])) the two anti-quarks, which are inclusive, can be identical as well as non- 
identical. Notice that in the above reactions the detected quark can be replaced by the detected 
anti-quark so that in reaction (|T7j ) one can also distinguish between the final states containing 
identical quarks and non-identical quarks. After mass factorization and renormalization for 
which we have chosen the MS-scheme one obtains the order a 2 contributions to the coefficient 
functions which are presented in ||. 

3 Notice that we make a distinction in nomenclature between the fragmentation densities and the 

fragmentation functions F k . 
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3 Review of the most important results 



The most important results of our calculations can be summarized as follows. From eq. ( fLOf) and 
the coefficient functions originating from the processes above we can obtain and ot corrected 
up to order a 2 



Me 2 ) = ^.'w 2 )[i + (^ 

+c A c F { 



'(C,{6} 



|,,6 C(3)-^} + n,C F T,{m(Z) + \} 
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(18) 



a L (Q 2 ) = ^ t (Q 2 ) 



ra s (/i 2 



47T 



[ C F {3} + 



a s (n 2 



47T 



M - y} 



+C A C F { 



11 In 



24 2023, r Q 2 74 



3» 



(19) 



Addition of <7l and or yields the well known answer for ctot © which is in agreement with the 
literature ||. Hence this quantity provides us with a check on our calculation of the longitudinal 
and transverse coefficient functions. Notice that in lowest order <r tot only receives a contribution 
from the transverse cross section whereas the order a s contribution can be only attributed to 
the longitudinal part. In order a 2 both o L and a? contribute to cr tot . 

Because of the high sensitivity of expression (^) to the value of a s , the longitudinal cross section 
provides us with an excellent tool to measure the running coupling constant. To illustrate the 
importance of the order a 2 contribution to <jl we have computed the ratio 



R L (Q 2 ) 



<tl(Q 2 



CTtotiQ 2 ) ' 

for fi = Q = M z and a s (5, M z ) = 0.126. The result is 

R L = 0.040 + 0.014 = 0.054 , (0.057 ± 0.005) 



(20) 



(21) 



where the first and the second number represent the order a s and order a 2 contribution respec- 
tively. Between the brackets we have quoted the result from OPAL 0. Here one observes a 
considerable improvement when the order a 2 contributions are included. Recently DELPHI 
used Eq. (ITUl) to determine the strong coupling constant. Their measurement of Rl yields 



R L = 0.051 ± 0.01 (stat.) ± 0.007 (syst.) , 
from which the strong coupling constant can be extracted. The result is 
af LO (5, M z ) = 0.120 ± 0.002 (stat.) ± 0.013 (syst.) . 



(22) 



(23) 



If one includes power corrections to ctl which are due to higher twist contributions of the order 
A/Q (see J7|) then one obtains 



a 



NLO+POW 



(5, M z ) = 0.101 ± 0.002 (stat.) ± 0.013 (syst.) ± 0.007 (scale) , 



(24) 
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where the scale uncertainty comes from varying the renormalization scale in he range 0.5 Q < fi < 
2 Q. The result above shows that the power corrections enhance 07, so that the strong coupling 
constant decreases. A similar analysis (see ||) shows that also the longitudinal fragmentation 
function Fl receives large order a 2 s contributions. However the corrections to the transverse and 
asymmetric fragmentation functions are small. 

Finally we discuss the order a 2 corrections to the asymmetric cross section. The latter is 
given by 



MQ 2 ) = E [ d ^ A 2 Qn = Qa [ ^<35(*. tfV) > ( 25 ) 



2^ rl 



Up to NNLO he above expression yields the following result 



a A = Q A [l + (^^) 2 { - 12/3oC(3)}] , A) = y Ca ~ t T f n f ■ 



3 3 

(26) 



Notice that a a is independent of the factorization scale. However it depends like <jl and a a on 
the renormalization scale via a s . From Eq. (|26| ) we infer that the NLO correction is zero but 
the NNLO contribution is nonvanishing although numerically it is very small for fi = Q = Mz- 
Finally in || we have evaluated the first and second moment of the asymmetric structure function 
and compared it with the data. The NNLO result for the first moment is 

f 1 dxF A (x, Mf ) = -0.016 , -0.0229 ± 0.0044 (OPAL) , 
J0.1 

-0.028 ± 0.006 (DELPHI) . 



(27) 



The second moment becomes 



C dx-F A (x, Ml) = -0.0020 , -0.00369 ± 0.00046 (OPAL) , 
J 0.1 2 

-0.0036 ± 0.0008 (DELPHI) . 



(2* 



We also computed the above sum rules up to NLO. However there is hardly any difference 
between NLO and NNLO which could already be expected from the comments made below 
Eq. (|26|) . The numbers presented above reveal that the experimental values are far below the 
theoretical predictions. Furthermore it turns out that the LO results are much closer to the ones 
obtained from experiment. In the case of Eq. ( p7|) we obtain —0.023 whereas for Eq. (p8|) we 
get —0.0027. Probably higher twist corrections might become very important as we saw in the 
case of &L below (E3|) . 
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